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Abstract  iai  i 

Bias  and  variance  equations  are  presented  for  two-dimensional  location 
estimators  of  a  nonmoving  point  source  of  radiation  in  an  isotropic,  sta¬ 
tionary  random  medium.  The  estimators  are  calculated  from  spatially  corre¬ 
lated  angle-of-arrlval  data  which  are  collected  simultaneously  at  two  sensor 
positions  and  assumed  to  consist  of  true  (unbiased)  source  angles  plus  zero- 
mean  angular  noise  with  equal  variances  at  both  sensors  and  negligible  higher 
moments.  Under  these  assumptions  the  square  of  the  estimator  bias  is,  in 
general,  a  quadratic  function  and  the  estimator  variance  a  linear  function  of 
the  spatial  data  correlation  coefficient.  However,  for  source  ranges  much 
larger  than  sensor  separation,  both  the  bias  and  the  variance  tend  to  increase 
linearly  with  decreasing  correlation  coefficient,  whereas  they  tend  to 
decrease  with  Increasing  sensor  separation.  The  combined  effect  for  a  distant 
source  in  a  stationary  random  medium,  when  evaluated  for  typical  spatial 
wavefront  autocorrelation  functions,  is  a  significant  reduction  in  the  esti¬ 
mator  bias  and  variance  dependence  on  sensor  separation,  as  compared  to  the 
uncorrelated  case.  With  minor  modifications,  the  same  results  apply  to  the 
equivalent  problem  of  using  time -of -arrival  data  from  three  colinear  sensor 
positions. 


Introduction 

In  sonar  or  radar,  passive  source  localization  is  concerned  with  the 
estimation  of  an  object's  location  from  its  emitted  random  radiation  (acoustic 
or  electromagnetic  waves).  For  the  two-dimensional  case,  this  estimation  of 
the  source  coordinates  (x,y)  requires  simultaneous  wavefront  angle-of-arrival 
measurements  by  at  least  two  separate  sensors  (or,  equivalently,  time-of- 
arrival  measurements  by  at  least  three  sensors).  In  many  cases  of  practical 
interest,  the  angle  measured  by  each  sensor  can  be  regarded  as  a  random  varia¬ 
ble  whose  mean  is  the  true  (unbiased)  source  direction  and  whose  variance  is  a 
constant  (independent  of  source  direction).  In  addition,  the  two  random 
variables  thus  generated  are  jointly  characterized  by  a  spatial  correlation 
coefficient  which  accounts  for  such  factors  as  wavefront  coherence,  sensor 
assembly  rigidity,  etc.  For  sound  propagating  through  the  ocean,  typical 
spatial  coherence  lengths  range  from  1  to  over  103  meters. M 

It  Is  well  known  that  even  though  the  angles  measured  by  the  sensors  may 
have  unbiased  means  and  constant  variances,  the  source  location  estimated  from 
these  angles  by  direct  calculation  has  a  bias  and  variance  which  are  functions 
of  the  source-sensor  geometry.  Several  authors  have  Investigated  this  problem 
under  the  assunption  of  spatially  uncorrelated  angle-of-arrival  (or  time-of- 
arrival)  data.  (2»3»4»5)  jn  paper.(  vie  present  equations  for  the  estimator 
bias  and  variance  resulting  from  spatially  correlated  data  (Section  I). 

Section  II  examines  the  bias  and  variance  dependences  on  correlation  coeffi¬ 
cient  and  source-sensor  geometry  and  presents  some  of  the  results  graphical¬ 
ly.  In  Section  III,  the  joint  effects  of  spatial  data  correlation  and  sensor 
separation  are  examined  for  distant  sources  In  a  stationary  random  medium. 
Conclusions  are  drawn  in  Section  IV. 


I.  General 


In  this  section,  we  develop  equations  for  the  bias  and  variance  of  source 
location  estimators  calculated  from  correlated  angle-of -arrival  data.  The  two 
sensors  are  positioned  at  x  *  ±D  In  the  x,y  plane  (Figure  1)  and  the  true 
source  location  Is  (xs,ys)  and  (R,8)  In  Cartesian  and  polar  coordinates, 
respectively.  The  true  source  angles  at  sensors  I  and  II  are  6}  and  62*  All 
angles  are  measured  counterclockwise  from  the  positive  x  direction. 

In  terms  of  0^,  82*  and  0  the  true  source  location  coordinates  are  as 
follows: 


R 


8 


( 


1+cos2(82-8 j^^cosfeg-e^Jcosfeg-^}  ^ 2 


tan 


sin  (02-e^) 
^  /  cos^-e^-cos^-^) 
V  sTnfepep 


) 


/  si n(02+O,)  \ 
xs  s  0  (  slnfeg-Sj)  / 

(costep-e.J-costep+fli) 
sin(82-01) 


(1) 

(2) 

(3) 

(4) 


(Note  that  these  and  the  following  equations  are  written  In  terms  of  (82+01) 
and  (02-0 1 )  to  simplify  their  Interpretation  for  specific  sensor-source  orien¬ 
tations.  Thus,  for  a  "broadside"  source  located  along  the  y-axls,  81+82  *  »•) 
Now  assume  that  the  angle-of -arrival  data  generated  at  sensor  positions  I 
and  II  are  random  variables  with  means  81  and  82  (i.e.,  true  source 

angles),  variance  o2  (same  for  both  sensors  and  independent  of  direction), 
correlation  coefficient  p  and  negligible  higher  moments.  The  source  coordi¬ 
nates  (x,y)  which  are  calculated  from  these  data  by  (3)  and  (4)  will  then  also 
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Figure  1.  Source-Sensor  Geoaetry. 


be  random  variables  whose  moments  can  be  determined  in  terms  of  the  angle-of- 
arrival  moments  by  expanding  (x,y)  into  a  Taylor  series  about  the  true  source 
location  (x$,ys).(6) 

Thus,  the  mean  and  variance  of  x  are  given  by 
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where  all  terms  involving  moments  hitter  than  o  and  p  have  been  neglected. 
Carrying  out  the  indicated  operations  on  (3)  we  obtain 
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Similarly,  the  mean  and  variance  of  y  are 
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\  * — r— 2 -  ([1  -  cos(e9+e,  Jcosfe^-e, )]2  +  sinz(e9+01)sin2(99-01] 

y  sin  (Sg-Bj)  \  c  1  4  1 

-  p[cOS(02-01)  -  (8) 


The  covariance  of  x  and  y  is  (again  retaining  only  terms  up  to  and  p) 
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Performing  the  Indicated  operations  on  (3)  and  (4)  yields 
2a2D2s1 0(02+0.)  / 

c  * - j -  (  COSfOg^)  -  005(02+0^0052(02-0^ 

sin  (Og-Sj)  \ 

+  pjcosfag^)  -  cos^-e^  (9) 

It  Is  convenient  to  consider  the  calculated  source  coordinates  (x,y)  as 
the  rectangular  components  of  a  two-dimensional  source-location  estimator 
(l.e. ,  a  two-dimensional  random  vector)  with  mean  (x,y)  and  variance  S.  From 
(5)  and  (7),  the  magnitude  B  of  the  estimator  bias  (i.e. ,  the  Euclidean  dis¬ 
tance  between  the  estimator  mean  (T,y)  and  the  true  source  location  (xs,ys)) 
is  obtained  as 

• 

B  -  ((x  -  Xj2  +  (7  -  yj2)1/2  -  a(e  +  PF  +  p2g)1/2  (10) 
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E  *  cos2(e2-e1)^l  -  2cos(0  2+fl1)cos(0  2-01)  +  cos2(e2-Q1)j 
F  *  2cos(e2-e1)^cos(e2461)(l  +  cos2(e2-e1))  -  2005(62-8^^ 
G  3  1  -  2cos(e2+e1)cos(82-e1)  +  cos2(92-e1) 

The  angle  6  of  the  estimator  bias  is,  from  (5)  and  (7): 
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Similarly,  from  (6)  and  (8),  the  estimator  variance  S  (i.e.,  the  mean- 
squared  Euclidean  distance  between  the  estimator  (x,y)  and  the  estimator  mean 
(7,y))  is  obtained  as 


S  3  a2  +  o2  3  H[I  +  pj] 


where 
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1 1.  Estimator  Bias  and  Variance 

The  strong  estimator  dependence  on  the  data  correlation  coefficient  p  Is 
shown  qualitatively  in  Figure  2  for  a  particular  source-sensor  geometry.  The 
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Figure  2.  Source-Location  Es'  imator  01*  rlbytlons  (shaded)  for  (a)  Positively 
Correlated  Data  (p  \),  //  Uncorrelated  Data  (p  ■  0),  and  (c) 

Negatively  Correlated  eta  (p  ■  -l). 
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cases  p  *+l  or  p  *  0  arise,  for  example,  when  the  angle-of -arrival  errors  are 
primarily  due  to  the  random  structure  of  the  propagation  medium  between  source 
and  sensors,  and  the  sensor  separation  is  much  smaller  or  much  larger,  respec¬ 
tively,  than  the  medium's  spatial  coherence  length.  Similarly,  the 
cases  p  *+l  orp  * -1  can  arise  when  the  errors  are  primarily  due  to  the 
random  motion  (rotational  or  translational)  of  the  rigid  two-sensor  assembly 
relative  to  the  propagation  medium. 

To  investigate  the  estimator  bias  and  variance  dependences  on  p  and  on 
the  source-sensor  geometry  In  more  detail,  we  make  use  of  the  inverse  rela¬ 
tions 


(12) 

(13) 


(Note  that  for  source  locations  In  the  first  quadrant,  both  02~®i  and  02+0l 
range  over  (0,ir)).  Substitution  of  (12)  and  (13)  into  (10)  and  (11)  yields, 
after  some  straightforward  but  tedious  algebra,  the  results  for  B,  6,  and  S: 
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Consider  first  the  estimator  bias  B.  Successive  differentiation  of  (14) 

M 

reveals  that  B  is  a  minimum  when  p  *  Pm-j n  =  -  "2n  *  or 


In  Figure  3,  pm^n  is  plotted  against  the  normalized  source  coordinates 

X  y  ^  y 

ClT  •  tt)*  where  (ff)2  *  (/)2  +  (b1)2  from  Fi 9ure  1-  We  note  that  Pm1n  >  0 

R2 

for  all  source  locations  except  those  in  the  region  cos2d  <  — *  <  1. 

D2 

By  substituting  (16)  Into  (14),  we  obtain  the  minimum  estimator  bias 

■w  ■  »»(i  -  - 
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B. 


mi  n 
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Sln26(l 


(17) 


where  B_  *  ^he  estimator  bias  obtained  for  uncorrelated  data,  i.e. 


for  p  *  0.  In  Figure  4, 


mi  n 


is  plotted  against  (jp  ,  -|). 


the  maximum  estimator  bias  BmflX  is  obtained  by  setting  p  *  +1  or  p  *  -1 


in  (14),  according  as  to  whether  pm^n  <  0  or  pm^n  >  0,  respectively.  The 
result  is 
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;  elsewhere 


(18) 
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X  y 

which  is  plotted  in  figure  5  against  (jp  ,  -jp) . 

Turning  now  to  the  estimator  variance  S,  we  note  from  (15)  that  S  is  a 

R2 

linearly  decreasing  function  of  p  when  — >  1  and  a  linearly  increasing  func- 
^2  D* 

tion  of  p  when  — 7  <  1.  Evaluating  (15)  for  p  *  -1  and  p  «  +1,  we  obtain  the 
u 

maximum  and  minimum  variances,  Smax  and  Sm^n: 
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where  SQ  *  P.Q  Is  the  estimator  variance  obtained  for  uncorrelated  data,  i.e, 


for  p  *  0.  Thus, 


S  S  . 

and  -Sr-E-  are  concentric  circles  In  the 


S™*Vr2 


(■»P  ,  •*=•)  plane.  Also,  or  ^  ,  according  to  whether  the  source 
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location  is  outside  or  Inside  the  circle  ^-k  •  1,  respectively. 

0£ 
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Finally,  It  should  be  noted  that  since  (14)  to  (19)  are  given  In  terms  of 
the  source-sensor  parameters  R,  0,  and  e (rather  than  the  measured  angles  6} 
and  d2),these  results  apply  equally  well  to  the  equivalent  problem  of  using 
spatially  correlated  tlme-of -arrival  measurements  from  three  collnear  sensors 
located  at  x  *  -D,  0,  and  +D  (Figure  1).  In  such  cases,  the  angular  noise 

O 

variance  a  used  here  must  be  evaluated  from  the  corresponding  time-delay 
variance  which  characterizes  the  difference  In  wavefront  arrival  times  between 
adjacent  sensor  pairs.  (An  example  Is  given  In  Section  III.) 


III.  Effect  of  Sensor  Separation  on  Distant  Source  Localization 

For  the  case  of  a  distant  source  (R  »  D)  we  obtain  from  (14)  and  (15) 
the  approximations 
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Inspection  of  (16)  to  (19)  (or  Figures  3  to  5)  shows  that  In  this  case 
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In  passing,  we  note  that  for  the  uncorrelated  case  (p  *  0),  (20)  and  (21) 

agree  with  the  corresponding  bias  and  variance  expressions  developed  In  refer 

2 

ence  4,  provided  that  the  time-delay  variance  oT  used  there  Is  converted  Into 
the  corresponding  angular  variance  used  here,  according  to 


2co 

0  "ffSfnl 


(where  c  ■  propagation  velocity). 

From  (20)  and  (21)  It  Is  seen  that  for  a  fixed  data  correlation  coeffi¬ 
cient  p  both  the  estimator  bias  and  variance  increase  with  decreasing  sensor 
separation  0,  whereas  for  a  fixed  sensor  separation  they  both  decrease  with 
Increasing  correlation  coefficient.  These  two  effects  therefore  counteract 
each  other  In  cases  where  p  is  a  decreasing  function  of  D.  Such  situations 
arise,  for  example,  when  the  angle-of -arrival  errors  are  primarily  due  to  the 
random  structure  of  the  propagation  medium  between  source  and  sensors. 

To  gain  some  Insight  Into  this  tradeoff  between  p  and  D,  we  represent 
both  (20)  and  (21)  by  the  functional  relation 


f(D,p)  »  (1  -  p)  (22) 

D 

where  K  Is  Independent  of  D  and  p.  The  net  effect  of  a  small  change  In  sensor 
separation  Is  therefore  given  by  the  derivative 
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dD 


(23) 


For  the  special  case  where  the  data  are  uncorrelated  for  all  sensor 
separations  we  have  p  •  0,  ^  »  0,  and 
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For  the  more  general  case  we  treat  the  angle-of-arrival  data  as  a  sta¬ 
tionary  random  process  with  spatial  autocorrelation  coefficient  p(D),  where 
p (0)  =  1,  p(»)  =  0,  and  is  restricted  by  the  requirement  that  p(D)  must  be 
positive-definite  (i.e.,  its  Fourier  transform  must  be  nonnegative  for  all 
spatial  wavelengths).  Since  this  restriction  cannot  be  expressed  in  a  conven¬ 
ient  algebraic  form,  we  evaluated  (23)  for  some  typical  spatial  autocor¬ 
relation  coefficients  and  different  values  of  p(D).  In  Table  I,  the  results 
are  normalized  relative  to  the  uncorrelated  case. (24).  The  fact  that 

^  <  0  for  all  cases  considered  suggests  that,  in  practice,  an  increase 
in  sensor  separation  will  always  result  In  a  net  decrease  of  estimator  bias 
and  variance  for  distant  targets.  However,  It  is  seen  from  Table  I  that  the 
rate  of  this  decrease  may  be  much  smaller  for  the  case  of  spatially  correlated 
data  than  for  uncorrelated  data. 

IV.  Conclusions 

We  have  Investigated  the  effect  of  spatial  correlation  between  two-sensor 
angle-of-arrival  data  on  the  source  location  estimates  calculated  from  these 
data.  The  conclusions  are: 

1.  Both  the  estimator  bias  and  variance  are  strong  functions  of  data 
correlation  and  source-sensor  geometry. 

2.  For  distant  sources,  estimator  bias  and  variance  both  tend  to  de¬ 
crease  linearly  with  Increasing  correlation  coefficient. 

3.  Data  correlation  dependence  on  sensor  separation,  as  given  by  typi¬ 
cal  spatial  autocorrelation  coefficients  in  stationary  random  media, 
can  cause  greatly  reduced  estimator  bias  and  variance  dependences  on 
sensor  separation. 


v:  ■ 
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4.  The  general  estimator  bias  and  variance  equations  (10,11,14,15)  are 
limited  in  accuracy  by  the  assumption  that  the  angle-of -arrival  data 
distributions  have  negligible  moments  above  the  second. 
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0.002 


0.005 


0.050 


0.053 


0. 


0.8 

0.013 

0.021 

0.10 

0.11 

0.7 

0.030 

0.050 

0.15 

0.18 

0.6 

0.055 

0.094 

0.20 

0.25 

0.5 

0.091 

0.15 

0.25 

0.33 

0.4 

0.14 

0.23 

0.30 

0.42 

0.3 

0.20 

0.34 

0.35 

0.52 

0.2 

0.27 

0.48 

0.40 

0.64 

0.1 

0.31 

0.67 

0.45 

0.78 

Table  I:  Estimator  Bias  and  Variance  Dependence  on  Sensor 


Separation  for  the  Following  Autocorrelation 
Coefficients  (a  >  0): 

a)  p(d>  .alapBl 

(2)  o(D)  -  exp(-aV) 


(3)  p(D) 


/  l-a|D|  ;  |D| 

1  0  ;  |D|  >i 
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